A possible determination of the quark radiation length in cold nuclear matter 
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We calculate the differential Drell-Yan production cross section in proton-nucleus collisions by 
including both next-to-leading order perturbative effects and effects of the nuclear medium. We 
demonstrate that dilepton production in fixed target experiments is an excellent tool to study 
initial-state parton energy loss in large nuclei and to accurately determine the stopping power of 
cold nuclear matter. We provide theoretical predictions for the attenuation of the Drell-Yan cross 
section at large values of Feynman xf and show that for low proton beam energies experimental 
measurements at Fermilab's E906 can clearly distinguish between nuclear shadowing and energy 
loss effects. If confirmed by data, our results may help determine the quark radiation length in cold 
nuclear matter Xq ~ 10~ 13 m. 
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I. INTRODUCTION 

The energy loss of a charged particle as it traverses 
dense matter is a fundamental probe of the matter prop- 
erties. Accurate theoretical calculations and experimen- 
tal measurements of this quantity became one of the great 
early successes of the classical and quantum theories of 
electromagnetic interactions 0. The stopping power, 
dE/dx, in the limit of large energies is dominated by 
bremsstrahlung processes and is related to the radiation 
length, Xq, of charged particles in matter as follows: 



dE/dx = E/X 



(1) 



In Eq. ([T]) Xq is approximately independent of the inci- 
dent particle momentum. Precise knowledge of the stop- 
ping power and radiation length of materials is widely 
used today in X-ray tomography, muon and proton ra- 
diography, and high energy nuclear and particle physics 
detector development and instrumentation. 

The fundamental constituents of nuclei, quarks and 
gluons, predominantly interact via their color charge. 
The forces between them are described from first princi- 
ples by the theory of strong interactions, Quantum Chro- 
modynamics (QCD), and are much larger in magnitude 
than the electromagnetic force. In the last decade, ad- 
vances in high energy many-body QCD have enabled ex- 
ploration of parton energy loss |2|-|5| in a novel hot and 
dense state of nuclear matter - the quark-gluon plasma 
(QGP). The predicted suppression of energetic particle 
production in nucleus-nucleus (A+A) collisions, dubbed 
jet quenching, is now definitively established [E0]- 



Alongside the excitement of this new discovery comes 
the realization that one of the biggest gaps in our current 
knowledge of nuclear reactions in extremis is the stopping 
power of cold nuclei for color-charged particles. In d+Au 
collisions at the Relativistic Heavy Ion Collider (RHIC), 
where the QGP is not formed, a similarly large suppres- 
sion of particle production at large forward rapidity is ob- 
served @, . STAR experimental data on ir° attenuation 
at y = 4, for example, can be fitted with disparate mod- 
els that emphasize either large shadowing or large nuclear 
stopping. A more realistic calculation that includes the 
Cronin effect (To| , hig h- twist shadowing and initial- 
state energy loss [l2[ - all independently constrained in 
different processes and center of mass energies - is pre- 
sented in Figure [TJ The incorporation of these effects in 
the perturbative QCD calculation is briefly summarized 
in Ref. [l^]. Our new results show that both coherent 
and inelastic scattering on the nucleus have comparable 
impact on the observed cross section attenuation in this 
kinematic range. This finding stresses once again the 
need to identify and reliably evaluate experimental ob- 
servables that are clean and sensitive signatures of jet 
energy loss in cold nuclear matter. 

From a phenomenological point of view, the sig- 
nificance of cold nuclear matter energy loss is well- 
established (T^ - fl7| . Theoretically, until recently, only 
final-state inelastic interactions in large nuclei relevant to 
semi-inclusive deeply inelastic scattering have been con- 
sidered [l8||. A new comparative study of initial-state 
and final-state energy loss [l2j suggests that in the limit 
of large parton energies these exhibit strikingly different 
path length and parent parton energy dependencies: 
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FIG. 1: Neutral pion suppression in minimum bias d+Au 
collisions at y = 4 at RHIC. Theoretical calculations that 
include known nuclear matter effects are shown. A complete 
simulation gives a good description of the experimental data. 



Here, Cr is the quadratic Casimir in the fundamental 
and adjoint representations for quarks and gluons, re- 
spectively, A s is the gluon mean free path of 0(1 fm) 
and a s is the strong coupling constant. In this paper 
we refer to the soft interactions with typical transverse 
momentum transfers squared £ 2 ~ 0.1 GeV 2 for quarks 
prior to the large Q 2 scattering as initial-state. Simi- 
larly, the soft interactions after the hard scattering are 
described as final-state. We emphasize that such sepa- 
ration is only possible if ^A 1 / 3 <C Q 2 . We will only 
be interested in lepton pair production of invariant mass 
squared M 2 = Q 2 > 10 GeV 2 , which is compatible with 
this constraint. Finally, we point out that the £ 2 de- 
pendence in Eq. © is implicit through k_lpm, the over- 
all suppression factor relative to the incoherent Bertsh- 
Gunion bremsstrahlung. 

The differential medium-induced bremsstrahlung spec- 
trum can be expressed as a solution of an inhomogeneous 
recurrence relation with suitably chosen boundary con- 
ditions [IH . These boundary conditions differ for initial- 
state and final-state energy loss and the results will be 
correspondingly different. Final-state interactions and 
Eq. (|3| have been investigated in detail, for example 
see [2l-Q- Let us now focus on Eq. ©. Our starting point 
is the integral form for the double differential medium- 
induced gluon bremsstrahlung spectrum (l2| : 
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In Eq. © k is the transverse momentum of the gluon 
relative to the direction of the parent parton, k + is 
its large lightcone momentum and q is the momentum 
transfer from the nuclear medium. The formation time 
of the gluon 77 = fc + /k 2 in comparison to the size of 
the medium L determines the degree of the destructive 
interference between the Bcrtsch-Gunion radiation and 



the radiation from the hard scattering. Let us focus on 
k + ~ E + and recognize that when tj L and k varies, 
the phase factor sin(L/r/) oscillates rapidly and aver- 
ages to zero. One is left with the first term in the inte- 
grand of Eq. ((4]), which is the incoherent medium- induced 
bremsstrahlung: 
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In the opposite limit ry 3> L for k + ~ E + we can expand 
the sine function to lowest order and obtain: 
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Note that the overall multiplicative coefficients {• • ■ } in 
Eqs. © and © have to be evaluated numerically with 
the relevant kinematic cuts specified in Ref. [l2|. In the 
coherent regime the coefficient also reflects the destruc- 
tive interference effect between the bremsstrahlung asso- 
ciated with the soft scattering and the bremsstrahlung 
associated with the large Q 2 process and can be numer- 
ically small. As the energy of the parent parton in the 
rest frame of the large nucleus grows, the approximation 
for A_E i r 1 ^: al _ stat0 given by Eq. © becomes more rele- 
vant. This is the basis for the advocated energy and 
path length dependence in Eq. ©. 

Of course, there are always parts of the emitted gluon 
phase space (fc + ,k) that are not compatible with sim- 
ple approximations. For this reason, we first evaluate 
the fully differential bremsstrahlung spectrum numeri- 
cally from Eq. ©, as described in [Tj|- From Eq. © in 
the small energy loss limit we can then quote a radiation 
length: 



Xn = LE 



dk A 



# kk+ dN°(E,L) 



dk+d 2 k 



(7) 



To summarize, for final-state interactions, the de- 
structive Landau-Pomeranchuk-Migdal (LPM) interfer- 
ence leads to a change in the functional form of radia- 
tive energy loss. Eq. © does not allow for a natural 
definition of a radiation length and implies that the ex- 
perimentally observable effects are limited to relatively 
small quark and gluon energies. In contrast, even if 
the LPM suppression factor Klpm ~ 1/10 in Eq. (J3J, 
AEl^ ial _ state retains some of the characteristics of inco- 
herent bremsstrahlung, see Eqs. ©, © and ©. For this 
reason, initial-state energy loss can also significantly af- 
fect experimental observables in he avy ion collider exper- 
iments of much higher ^snn [H, EH Hil . Furthermore, 
Eq. © implies that the stopping power of cold nuclear 
matter for partons prior to a hard Q 2 3> ^qcd scatter- 
ing can be characterized by a radiation length Xq defined 
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in Eq. ((7]) . One can see parametrically from Eq. <j2j) that 
X is expected to be of O(10 fm — 100 fm) - the short- 
est radiation length in nature, ten orders of magnitude 
smaller than the radiation length of high-Z materials, 
such as tungsten, for electrons. 

The Drell-Yan process in heavy ion collisions [I?} - 
q + q — > 7* — > / + + l~ at lowest order (LO) - is an ideal 
probe of initial-state effects. The final-state particles do 
not interact strongly with the nuclear medium, providing 
a relatively clean experimental signature. Still, definitive 
separation of leading-twist shadowing effects [2l| - |23l ] and 
parton energy loss fl2 , [HHItI has so far proven challeng- 
ing 0, [2t| . In no small part this difficulty arises from 
the fact that the very same Drell-Yan data in proton- 
nucleus (p+A) reactions is used to constrain the shad- 
owing parameterizations [30| . Going to higher center of 
mass energies, such as what is available at RHIC and the 
Large Hadron Collider (LHC), will only amplify this un- 
fortunate ambiguity. In this paper we demonstrate how 
dilepton production in conjunction with a suitably low 
proton beam energy -Ebcam = 120 GeV, corresponding to 
Y / s]v7V = 15 GeV, at Fermilab's experiment E906 can dis- 
entangle initial-state energy loss and nuclear shadowing 
effects and pinpoint the stopping power of cold nuclear 
matter for quarks and gluons with ~ 20% accuracy. 

Our work is organized as follows. In Section [TT] we 
investigate lepton pair production in hadron-hadron re- 
actions at next-to-leading order (NLO) in collinear fac- 
torized perturbative QCD and validate the theoretical 
simulation tools against experimental data. In Section lTOl 
we discuss the cold nuclear matter (CNM) effects: initial- 
state energy loss and nuclear shadowing. A complete cal- 
culation of the dilepton production rate in p+A collisions 
at NLO that includes CNM effects is given in Section ITVl 
We compare results to existing data on the attenuation 
of the Drell-Yan cross section in reactions with nuclei 
and present predictions for this suppression versus the 
cold nuclear matter radiation length Xq for the upcom- 
ing E906 measurements. A summary and conclusions are 
presented in Section IVl 



II. THE DRELL-YAN PROCESS IN P+P 
COLLISIONS 

In this paper we will be interested in the differential 
Drell-Yan production cross sections da DY jdxpdQ 2 and 
dcr DY I dydQ 2 where Q 2 = q 2 = (p t + +Pi- ) 2 is the invari- 
ant mass squared, xf = ^qhj \fs is the Feynman xf, and 
V = \ ln[(<7o + 9l)/(<Zo ~ <7l)] is the rapidity of the lepton 
pair. For illustration, we start the discussion with the LO 
cross section in the collinear factorization approach (3l| : 



do 



DY 



E 



Aire^a 2 



dxpdQ 2 9Q 2 s(xi + X2) 



fq(q){xi,li)fq(q)(x2,fi) -(8) 



q,q 



In Eq. ©aw 1/137, e q is the quark fractional electric 
charge and fi(xi, (i) are the parton distribution functions. 




0.1 0.3 0.5 0.7 



4 6 8 10 12 14 
Q [GeV] 



FIG. 2: Next-to-leading order calculations for the muon 
pair production cross section in p+D collisions at y^ijvjv = 
38.8 GeV is compared to Fermilab E772 experimental data. 
The yellow bands represent the theoretical uncertainty due to 
choice of factorization and renormalization scales. We have 
only included isospin effects for the deuterium target. 



The incident parton lightcone momentum fractions x\ , x-i 
obey the relations: 



x\Xi = Q 2 /s , xf = xi - x 2 



(9) 



It is well-known that the perturbative next-to-leading 
order corrections to the Drell-Yan process can be quite 
significant [32T - l3~4l | . The NLO calculation also adds 
a Compton scattering contribution to the annihilation 
graphs. Schematically, the one- loop cross section can be 
written as: 



da DY 

NLO 



dog Y 



dQ 2 dxF dQ 2 dxF dQ 2 dxF 



(10) 



The importance of going to 0(a 2 a s ) is two- fold. First, 
it removes phenomcnological Knlo factors and allows a 
more reliable evaluation of the perturbative uncertainty 
through the variation of the renormalization and factor- 
ization scales fiR — (if = /«. Second, it facilitates a more 
accurate incorporation of dynamical nuclear effects that 
scale with the quadratic Casimir, which plays the role 
of the average squared color charge for quarks and glu- 
ons. In this work we implement the NLO Drell-Yan cross 
section calculation following Ref. (33|. 

An alternative approach to lepton pair production is 
based on the dipole model [35|, [36[. In this model the 
Drell-Yan process is viewed as induced virtual photon 
bremsstrahlung in the target rest frame and nuclear ef- 
fects are parameterized in the dipole cross section. In 
this paper we adhere to the traditional QCD factorization 
approach where higher-order corrections can be system- 
atically evaluated order-by-order in perturbation theory. 
Furthermore, nuclear effects in perturbative QCD are in- 
corporated at the cross section level as opposed to the 
amplitude level and are thus dependent on the average 
path length of partons through dense matter (48j . 
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To test the baseline perturbative calculation for the 
Drell-Yan production cross section we compare our simu- 
lation results to Fermilab E772 data [13] from a fixed deu- 
terium target experiment with an incident proton beam 
energy I?beam = 800 GeV in Figure |2l We take into ac- 
count isospin effects as described in [381 ] . The left panels 
of Figure [5] show results for two different values of the 
lepton pair mass Q. The right panels show results for 
two different values of xp. The yellow band represents 
the perturbative uncertainty due to the variation of the 
renormalization and factorization scales Q/2 < p < 2Q. 
We have checked that shadowing or energy loss effects for 
the deuterium target have less than a few % effect on the 
calculated cross section and are not shown in Figure [5] 
for simplicity. 



III. COLD NUCLEAR MATTER EFFECTS 

As seen in Figure [TJ cold nuclear matter effects can 
significantly alter the production cross section for ener- 
getic/massive final states. Universal initial-state leading- 
twist shadowing effects are incorporated in nuclear par- 
ton distributions (nPFDs). Application for different ob- 
scrvables or different v / stvjv is achieved by multiplying 
the standard parton distribution functions (PDFs) with 
a nPDF correction factor [38[ . Alternatively, one can at- 
tempt to calculate CNM effects associated with the elas- 
tic, inelastic and coherent initial-state and final-state par- 
ton interactions in large nuclei from first principles [4l|. 
This latter approach aims to elucidate the physics that 
underlays the observed nuclear modification. To date, 
there is no global nPDF analysis that attempts to sepa- 
rate the process-dependent and process-independent ef- 
fects, largely due to the limited knowledge of the latter. 
Better theoretical and experimental control on cold nu- 
clear matter energy loss is one step to help rectify this 
deficiency. 

The process-dependent nuclear effects in question are 
enumerated below. First comes the Cronin effect that 
is often modeled through initial-state transverse momen- 
tum broadening [42[ and affects the qp distributions of 
energetic particles. In this manuscript wc consider qp- 
integrated cross sections, see for example Eqs. © and 
(fTU|) . and do not include the Cronin effect. Next comes 
dynamical shadowing that arises from the coherent final- 
state interactions of the recoil parton in the nuclear tar- 
get [ll[ . For the Drell-Yan process at LO this effect van- 
ishes since there is no parton in the final state. At NLO 
there is indeed a parton in the final state, however, for 
the dilepton masses of interest Q = 4—10 GeV the power 
suppressed high-twist shadowing effects, cx (A 1 / 3 £ i 2 )/Q 2 
where £ 2 ~ 0.1 GeV 2 , are expected to be negligible. This 
leaves the stopping power of large nuclei for incoming 
quarks and gluons as the most significant dynamical nu- 
clear effect for Drell-Yan production. We discuss the cold 
nuclear matter energy loss below since it is central to our 
paper. 



We evaluate initial-state energy loss numerically for 
minimum bias p+A reactions with deuterium ( 2 D), 
beryllium ( 9 Be), carbon ( 12 C), aluminum ( 27 A1), iron 
( 56 Fe), and tungsten ( 184 W) targets following Ref. [ilj . 
We use A g = 1 fm for the gluon mean free path and vary 
the momentum transfer per interaction with the medium 
from £ = 0.175 GeV to £ = 0.5 GeV H|. Comparing 
our results to Eq. ([7]) we find that in the limit of small 
energy loss A^ i ' I ^ ial _ statc < E or, cquivalently, L <C X , 
the radiation length of 100 GeV quarks in cold nuclear 
matter is between 30 fm and 160 fm for the parameters 
discussed above. We also calculate the probability distri- 
bution P q g (e) for quarks and gluons to lose a fraction of 
their energy e = ^\ uii/E due to multiple gluon emission 
in the Poisson approximation (38l . |39| : 



f P g>g (e)de = 1, 
Jo 



(11) 



A TT-rad. 

7-j / \ , q.q initial — state ~\ 

P q , g (e)ede = ^_ . (i 2 ) 



The calculation and/or implementation of initial-state 
or final-state state energy loss processes focuses on real 
medium-induced gluon emission. Let us consider for 
definitiveness the energy loss of the first incoming parton 
(characterized by lightcone momentum pf and fraction 
xi). The differential distribution of emitted gluons itself 
is evaluated by separating it from the partonic scattering 
cross section. For a general hard process: 



da 
dPS 



dxidx 2 f{x\)f{x2) J de 
dN9(e) 



(|M[p+(l-e),p+ 
2.ti(1 — e)x2S 



+ 1|2\ 



de 



(13) 



In Eq. (|13p e = k + /pf and parton flavor indices and the 
sum over flavors are not shown explicitly. We emphasize 
that the energy loss results only depend on the details 
of the final state through the value of p\ and in the 5- 
function pi = pi(l — e) +p 2 - Wc have also integrated over 
the gluon's transverse momentum k, noting that the di- 
rection of the incoming parton does not change on aver- 
age. In the soft gluon approximation, which we use in this 
study, medium-induced bremsstrahlung can be identified 
even at the amplitude level 0, [l2j ■ Proper factorization 
of the differential energy loss distribution in the e — > 1 
limit can only be shown in QCD at the cross section 
level . This distribution does include the interference 
of the bremsstrahlung from the soft in- medium scattering 
with the bremsstrahlung from the large Q 2 process. 

For a single (n g = 1) emitted gluon dN 9 (e)/de = P(e) 
is also the probability distribution for fractional energy 
loss e. In general, n g ^ 1 and the probability distri- 
butions Pq t g{e), see Eqs. (fTTj) and (fl"2|) . is constricted 
from dN 9 (e)/de in the independent Poisson approxima- 
tion [38l . |39| to gluon emission. Changing variables, 
x\ — > xi/(l — e), and the order of the integration we 



find: 



da 
dPS 



deP(e)f 



1 



2xiX2S 



1 - e y 

(|A/[p+ p+]| 2 >(2vr)^%- P/ ).(14) 



From Eq. (fT4")l . initiai-state energy loss is most conve- 
niently implemented in the perturbative calculations of 
the differential Drcll-Yan cross section Eqs. © and (JTUJ) 
as follows: 



(15) 
(16) 



Here, /(a:, /z) is the distribution function of the par- 
ton from the incident proton. The physical meaning of 
Eqs. (fT5")) and (|TB| is that in the presence of nuclear stop- 
ping larger quark and gluon energies in the nuclear wave 
function are probed by the same lepton pair kinematics. 
These results imply that the attenuation of the Drell- 
Yan cross section depends not only on the magnitude of 
cold nuclear matter energy loss but also on the steepness 
of the PDFs. Specifically, large suppression is expected 
at large values of Feynman x p (or large rapidity) where 
f'(x, fi) <C 0. We have denoted by f'(x, fi) the derivative 
of the parton distribution function with respect to the 
momentum fraction x. Finally, we note that accounting 
for fluctuations in the energy loss via Pq, g {e) is critical 
when jet and particle production is near kinematic lim- 
its. We have checked that using the mean energy loss 
A ^™g d initiai-statc/^ instead of the full convolution over 
Pq t g(e), overpredicts the suppression at large xp by more 
than a factor of two. 



IV. LEPTON PAIR PRODUCTION IN P+A 
COLLISIONS 

We now present results from the full calculation of 
Drell-Yan production at NLO that include cold nuclear 
matter effects. We first consider existing measurements 
from the Fermilab E772/E886 experiment [l6| which 
were performed using several nuclear targets in the path 
of a proton beam of energy E\ a \> = 800 GcV. Dilepton 
data and theoretical simulations are presented in Fig- 
ure [3] as the ratio of the differential cross sections on two 
different targets in minimum bias p+A reactions scaled 
down to a binary nucleon-nucleon interaction: 



Ratio(A/B) = R 



DY 
AB 



B 



A- 



dxpdQ 2 ' dxpdQ 2 



(17) 



We considered two different invariant mass ranges 
4 GeV < Q < 5 GeV (left panels) and 6 GcV < Q < 
7 GeV (right panels). We also included two different 
ratios: JJ^, and in the top and bottom panels, 
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FIG. 3: Plot of the ratio of the dimuon cross section 
cr(p+W)/(T(p+A), where A = D,Be, from both experiment 
and theory for 800 GeV protons colliding with a fixed nuclear 
target. The theory curves are the NLO Drell-Yan cross sec- 
tion calculation, including shadowing effects and initial-state 
energy loss, respectively. The yellow bands indicate a five-fold 
variation of the mean quark medium-induced energy loss. 



respectively. The dashed lines represent a calculation of 
the cross section ratios that incorporates a leading twist- 
shadowing parameterization [3(J. The solid lines show 
our simulation with initial-state energy loss [l2[ that is 
compatible with the theoretical description of RHIC d+A 
and A+A data [U, M, HI ■ 

It should be noted that while both calculations de- 
scribe the E772 /E866 data with equal acumen, they em- 
phasize completely different physics. The former simula- 
tion assumes a modification of the nuclear wave-function 
and is driven by the decreasing momentum fraction X2 of 
the parton in the nucleus when Feynman xp grows. The 
behavior of R^b for the latter case is governed by the en- 
ergy loss of the incident parton that carries an increasing 
fraction x\ of the momentum of the proton. The current 
accuracy of the experimental data does not allow one 
to distinguish between these two competing paradigms. 
Furthermore, it cannot precisely constrain the magnitude 
of Ai5™^j al _ statc . The yellow band in Figure [3] represents 
a variation of the typical momentum transfers between 
the incoming parton and the medium £ in the range of 
0.175 GeV to 0.5 GeV, as described in Section EH Such 
variation induces an approximately five-fold uncertainty 
in the magnitude of the stopping power of large nuclei for 
incident partons. While the extreme cases of very large 
or very small energy loss are not favored by the data, 
a considerable uncertainty can still exist for this physics 
scenario. 

A unique opportunity to disentangle initial-state en- 
ergy loss effects from nuclear shadowing is presented by 
the approved and now under construction Fermilab ex- 
periment E906 [M EH] . It is a fixed target Drell-Yan 
experiment with -Bbeam = 120 GeV. Its primary goal 
is to accurately measure the large Xp d/u asymmetry 
and it can comfortably cover the 0.4 < xp < 0.9 and 



Mj/ip < Q < Mr kinematic range. A secondary goal 
of E906 is to precisely determine cold nuclear matter 
energy loss. Its low center of mass energy per nucleon 
pair y/SNN = 15 GeV allows a selection of the dilepton 
mass which practically eliminates shadowing effects. In 
contract, nuclear stopping is expected to be even more 
important at these energies. 



Next, let us review again the prerequisites for the de- 
termination of the stopping power of cold nuclear matter 
and its radiation length Xq. First, one needs to iden- 
tify the part of phase space where shadowing effects arc 
minimal (we are interested specifically in the Drell-Yan 
process). Second, experimental measurements must es- 
tablish a statistically significant suppression of the dilep- 
ton production cross section in p+A reactions R9][ < 1 
at large xf- Third, one has to confirm an approximately 
linear dependence of R®X on the nuclear size A 1 / 3 . Note 

that the A 1 / 3, dependence of the AE lad - is correlated with 
a linear dependence on the parton energy. This will allow 
for a proper extraction of X , 
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FIG. 4: Next-to-leading order theoretical predictions for the 
Drell-Yan dimuon cross section attenuation cr(p+A)/cr(p+D) 
at Fermilab's E906. As with Figure[3] the curves include shad- 
owing effects (dashed line) or initial-state energy loss (solid 
lines), respectively. We have considered three different radia- 
tion lengths Xq of quarks in large nuclei, ranging from 30 fm 
to 160 fm. The lower right panel highlights the difference be- 
tween the two physics scenarios at large xf as a function of 
A 1/3 . 



Our numerical results, relevant to the E906 program, 
are presented in Figure |H We have chosen to consider 
da DY / dxpdQ 2 versus Feynman xp for Q = 4.5 GeV, 
and have carried out simulations for deuterium, carbon, 
iron, and tungsten targets. In order to obtain the ratio 
defined in Eq. (|17[) , we have normalized the nuclear cross 
sections to that of a D target which minimizes any trivial 
isospin dependence. In all cases, the effects of shadow- 
ing, which are shown with dashed lines in Figure HI are 
negligible. On the other hand, energy loss effects, repre- 
sented by solid lines for three different radiation lengths 
Xq = 30 fm, 50 fm and 160 fm, can clearly be detected, 
especially for large values of xf- We have also included 
in Figure [4] simulated experimental data which shows the 
anticipated E906 statistical precision for R^X D ■ With 
careful selection of the heavy nuclear targets, such as Fc 
and W, and the lepton pair kinematics, the two com- 
peting pictures for the physics that underlays the sup- 
pression of dilepton production in p+A reactions can be 
definitively confirmed or refuted. As seen in Figure |4j 
the radiation length of cold nuclear matter for quarks 
can likely be constrained with ~ 20% accuracy. 

The bottom right panel in Figure |4] shows the atten- 
uation of the Drell-Yan cross section at a fixed value of 
Xf = 0.9 versus the linear size ~ A 1//3 of the target nu- 
cleus. Our results reflect the fact that initial-state energy 
loss of quarks and gluons prior to a large Q 2 scatter- 
ing depends approximately linearly on the path length 
through strongly interacting matter. This is in contrast 
to final-state energy loss where quadratic path-length de- 
pendence 0-13] would suggest A 2 / 3 scaling of the attenu- 
ation of the lepton pair production rate. E906 measure- 
ments will also help clarify this open question. 



V. SUMMARY AND CONCLUSIONS 

One of the most pressing current questions in the the- 
ory and phenomenology of heavy ion reactions at high 
energies is related to the stopping power of cold nu- 
clear matter for energetic quarks and gluons. Plenti- 
ful experimental evidence and theoretical justification 
exist for the significance of initial-state energy loss ef- 
fects [I, H GJ, III El- It is, therefore, unfortu- 
nate that attempts to precisely determine their strength 
have so far been inconclusive. These earlier studies em- 
ployed a functional form for the energy loss that might 
have been better-suited to describe final-state scattering 
processes in the QGP, and yielded large uncertainties for 
-dE/dx [H, |H [H, H|]. Little is known quantitatively 
for the interactions of partons and hadrons in cold nu- 
clear matter beyond the interaction lengths, or mean free 
paths, for selected few species, such as protons or al- 
pha particles. Only recently have developments in the 
theory of quark and gluon propagation in cold nuclear 
matter suggested that initial-state energy loss in QCD 
retains certain characteristics of induced electromagnetic 
bremsstrahlung [l2j |. Thus, it can be described in terms 
of a radiation length Xq. A real opportunity exists today 
to determine this shortest radiation length in nature. 

To this end, we have embarked on a program to de- 
veloped the theoretical tools that can both motivate 
and facilitate the upcoming E906 experimental p+A pro- 
gram. To ensure that the baseline differential Drell-Yan 
cross sections in elementary nucleon-nucleon reactions 
are reliably estimated, we opted for a next-to-leading or- 
der accuracy of the perturbative QCD calculation. The 
new simulation tool was validated against Fermilab's 
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E772 muon pair production measurements in p+D reac- 
tions. Next, we implemented initial state-nuclear effects, 
such as leading-twist shadowing and medium-induced 
bremsstrahlung in large nuclei. Care was taken to in- 
corporate the fluctuations in the magnitude of parton 
energy loss. We showed that neglect of these fluctuations 
will result in a sever overestimate of the attenuation of 
the lepton pair production rate. Results from our work 
have been used to evaluate the dilcpton background for 
the T measurements at RHIC via the l + + l~ decay chan- 
nel m. 

For nuclear targets with large mass number we 
first verified that existing Drell-Yan measurements can- 
not distinguish between calculations that incorporate 
initial-state energy loss and leading-twist shadowing 
paramctrizations. While this observation has been made 
before, our work is more general in that it goes beyond 
the fixed energy loss per unit length conjecture. Our 
results also suggest that moving to the higher center-of- 
mass energies of RHIC and the LHC will not help sepa- 
rate the two competing physics effects. We showed that 
radiation lengths in cold nuclear matter for quarks Xq 
of the order of 50 fm to 100 fm are compatible with the 
£bcam = 800 GeV Fcrmilab E772/E866 data. The same 
stopping power of large nuclei has worked well for the the- 
oretical description of hard probes production at RHIC. 

To definitively confirm or refute initial-state energy 
loss as the principal nuclear effect which leads to the 
suppression of the Drell-Yan cross section in p+A reac- 
tions we presented predictions for the upcoming Fcrmilab 
experiment E906. We exploited its low center-of-mass 
energy to identify a kinematic region where shadowing 



effects vanish. We demonstrated that, with the antic- 
ipated statistical accuracy of the data in the large xf 
region, radiation lengths as large as 200 fm can be de- 
tected. If Xo is of O(50 fm), there is a good chance that 
it can be determined with ~ 20% accuracy. 

As a final note, we recognize that experimental 
measurements may differ from the most sophisticated 
theoretical predictions. With the next-to-leading order 
simulation tools at hand, upcoming E906 experimental 
data can be quickly and reliably analyzed for parton 
energy loss effects. The theoretical progress described 
here, combined with experimental advances in the use of 
dimuons, is expected to result in the first una mbig uous 
measurement of quark energy loss in nuclei [46j and 
provide the much-needed standard candle for gauging 
the nuclear response to strongly interacting particles. 
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In this paper we will use the term heavy ion collisions 
to describe both p+A and A+A hadronic reactions at 
relativistic energies 

For example, in a picture of interfereing virtual pho- 
ton production amplitudes from distinct nucleon-nucleon 
scatterings the path length will be the average of the two 
path lengths. 

The momentum transfer between the jet and the medium 
is often denoted by \i. We use £ in this paper to indicate 
that the default value corresponds to the same dimen- 
sional scale used in dynamical shadowing calculations. 



